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摘要 : 利用 LAMOST (Large sky Area Multi-Object fiber Spectroscopic Telescope) DR7 和 


部 分 DR8 的 


分 辨 率 数据 与 Chen 等 人 四 的 WISE 食 双星 的 星 表 进 行 交 又 ， 得 到 了 151 颗 短 周 


期 的 太阳 型 恒 


星 ， 再 选择 其 中 观测 次 数 在 9 次 以 上 ， 并 且 相 位 分 布 比较 均匀 的 23 颗 星 ， 通 过 对 


分 析 ， 发 现 所 


它们 视 向 速度 


行 分 析 ， 得 到 了 这 批 样本 的 轨道 周期 、 质 量 比 、 轨 道 偏心 率 等 轨道 参数 。 通 过 
得 到 的 周期 与 WISE 给 出 的 周期 基本 一 致 ， 所 以 WISE 的 周期 是 可 靠 的。 之 后 / 


WISE 的 周 基 


作为 已 知 量 ， 再 去 求解 其 他 轨道 参数 ， 这 样 使 其 他 的 轨道 参数 更 加 准确 。 ARI 


拟 合 较 好 的 23 颗 食 双星 ， 对 这 批 参数 的 统计 性 质 进 行 分 析 ， 并 与 Ragavan 等 人 ™ 


其 中 速度 曲线 
的 结果 进行 对 


比 ， 发 现 与 Ragavan 等 人 "的 结果 基本 一 致 ， 在 周期 10 d 以 内 的 星 ， 质 量 比分 布 


在 0.2~1 范围 


内 P KLH M, AI Ah 心率 都 比较 小 ， 在 0.2 以 下 。 
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双星 在 各 种 


学 家 甚至 认为 绝 


的 分 类 方法 ， 根 
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不 同 质量 的 恒星 系统 中 都 很 常见 ， 从 星团 到 整个 星系 都 是 如 此 ， 一 些 天 文 
大 多 数 的 恒星 都 处 在 双星 系统 中 。 双 星 有 很 多 种 类 型 ， 由 Kopal" 介绍 
据 其 中 的 一 颗 或 两 颗 子 星 是 否 充满 洛 希 办 ， 一 般 可 以 分 为 分 离 双 星 、 半 相 接 


。 当 两 颗 星 均 未 充满 洛 希 辨 时 ， 两 颗 星 独立 演化 ， 这 被 认为 是 分 离 双 星 ;， 有 


双星 和 相 接 双星 
一 颗 星 充满 洛 希 瓣 是 半 相 接 双星 ;， 当 两 颗 星 都 充满 洛 希 瓣 时 是 相 接 双 星 ， 相 接 双星 两 颗 星 


= 


成 公共 包 层 并 发 生物 质 转移 ， 两 颗 星 相互 作用 很 强 。 双 星 的 质量 比 、 各 成 员 星 质量 、 轨 
道 等 其 他 参数 可 由 食 双 星系 统 光 变 曲线 的 分 析 以 及 时 域 光 谱 观 测 来 测定 。17 世纪 ， 现 代 天 
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文学 诞生 后 不 久 ， 人 们 就 认识 到 双星 系统 和 多 星系 统 的 存在 ，1935 年 Kuiper ”最 早 提 出 ， 
从 理论 上 看 ， 确 定 恒星 的 多 重 性 以 及 轨道 参数 的 分 布 非常 有 价值 。 恒 星 的 多 样 性 是 恒星 形成 
的 常见 结果 。 因 此 ， 研 究 双 星 比 和 双星 系统 的 主要 性 质 以 及 它们 对 初始 质量 和 环境 的 依赖 性 
是 检测 该 形成 过 程 的 有 力 工 具 。 恒 星 中 双星 比 的 表征 和 轨道 参数 的 统计 研究 对 于 天 体 物 理 
学 的 许多 领域 至 关 重 要 。 双 星 的 演化 不 仅 对 星 族 的 演化 有 很 大 影响 ， 而 且 可 以 解释 许多 重要 
的 天 体 物理 学 现象 ， 例 如 Ia 型 超新星 ， 两 个 黑洞 并 合 的 引力 波 事件 等 。 无 论 是 星团 还 是 整 
个 星系 ， 双 星 在 各 种 规模 的 恒星 系统 中 都 很 丰富 。 对 于 双星 的 统计 分 析 ， 有 利于 我 们 理解 恒 
星 在 不 同 环境 中 的 形成 与 演化 。 双 星 的 统计 性 质 可 以 由 双星 比 、 质 量 比分 布 、 轨 道 周期 分 布 
和 轨道 偏心 率 分 布 很 好 地 描述 出 来 ， 并 且 通 常 这 些 统计 性 质 还 与 恒星 群 的 年 龄 、 金 属 丰 度 以 
及 主星 质量 有 关 ” TS, 
研究 双星 有 很 多 种 方法 ， 比 如 用 光 变 曲线 的 方法 ， 用 视 向 速度 的 方法 ， 或 者 用 赫 罗 图 
的 方法 。 在 双星 系统 中 ， 由 于 它们 围绕 共同 质心 做 轨道 运动 ， 它 们 的 投影 径 向 速度 的 变化 ， 
可 以 在 轨道 周期 确定 后 ， 根 据 轨道 相位 绘制 出 来 。 所 以 对 这 样 的 速度 曲线 进行 分 析 ， 就 可 
= 以 确定 描述 轨道 几何 形状 的 量 以 及 与 两 颗 恒 星 的 质量 有 关 的 量 。 之 后 在 1890 4F Algol 轨道 
= 被 发 现 后 不 入， 就 在 此 基础 上 发 展 出 很 多 图 形 化 的 方法 ， 随 着 计算 机 的 出 现 和 发 展 ， 这 些 
方法 都 得 到 了 进一步 的 完善 和 改进 。 第 一 条 被 成 功 分 析 的 食 双 星 的 光 变 曲线 是 圆 形 轨道 的 
恒星 给 出 ， 当 时 的 主要 技术 是 由 Russell 和 Shapley 开发 出 的 方法 ， 后 来 经 过 Kopal 的 广泛 
使 用 并 进一步 发 展 ， 并 且 定 义 出 了 著名 的 双星 分 类 方案 。 这 之 后 慢 慢 探索 出 用 侍 里 叶 方 法 
分 析 光 变 曲 线 ， 尤 其 是 对 于 明显 的 非 球 形 的 恒星 系统 ， 同 时 还 受到 其 他 相互 作用 影响 的 系 
统 。 另 一 方面 ， 计 算 机 的 日 益 强 大 使 研究 人 员 能 够 考虑 在 表面 亮度 分 布 不 均匀 的 非 球 形 恒星 
围绕 其 共同 的 质量 中 心 旋转 时 产生 光 变 曲线 的 问题 ， 发 展 成 了 可 以 考虑 到 所 有 当时 已 知 天 
本物 理学 的 数值 模型 中 合成 双星 光 变 曲线 的 方法 。Lucy™””、Mochnacki 和 Doughty ™ 7, 
Rucinski 关于 相 接 双星 的 工作 ，Hill 和 Hutchings”, Hutchings 和 Hil 以 及 Wilson 
和 Devinney 关于 所 有 类 型 双星 的 工作 导致 了 对 于 双星 的 分 析 方 法 的 巨大 变化 。 他 们 编写 
出 很 多 解 双星 轨道 参数 的 程序 ， 这 些 程 序 极 大 地 提高 了 我 们 可 以 确定 多 种 类 型 的 双星 系统 
的 物理 特性 的 精度 ， 并 促进 了 对 一 般 双 星 演 化 模型 的 更 严格 的 测试 。 
近年 来 ， 随 着 大 规模 天 文 巡 天 观测 的 出 现 ， 我 们 获得 了 大 量 恒星 的 高 质量 测 光 、 光 谱 巡 
= 天 的 数据 。 如 LAMOST™, SDSs™™, 2MASS™, Gaia ”巡天 等 ， 对 于 这 些 大 样本 数据 ， 
= 发 展 出 了 针对 大 样本 数据 的 统计 方法 来 研究 双星 的 性 质 。 使 用 统计 方法 ， 不 需要 区 分 每 个 星 
是 单 星 还 是 双星 ， 只 需要 分 析 样 本 的 总 体 统计 性 质 。 因 此 ， 这 些 统计 方法 可 以 大 大 增加 样本 
量 ， 研 究 不 同样 本 的 双星 性 质 ， 如 金属 丰 度 、 有 效 温 度 、 双 星 比 、 轨 道 参 数 等 变化 。 
金属 丰 度 、 年 龄 、 质 量 等 都 与 恒星 的 形成 和 演化 密切 相关 ， 而 近年 来 很 多 工作 都 表明 这 
些 都 与 双星 的 统计 性 质 有 非常 密切 的 关联 。 例 如 ， 现 在 已 经 有 很 多 证 据 表 明太 阳 型 恒星 的 
密 近 双星 的 双星 比 随 着 金属 丰 度 的 增 大 而 明显 减 小 。Moe 等 人 "为 了 能 够 准确 地 测量 密 近 
双星 的 双星 比 与 金属 丰 度 的 关系 ， 将 一 些 有 分 歧 或 者 矛盾 的 结果 进行 了 整体 的 分 析 和 综合 ， 
在 校正 完备 性 后 ， 所 有 不 同 的 密 近 双星 数据 集 在 双星 比 和 人 金属 丰 度 之 间 都 表现 出 一 致 的 强 
反 相 关 性 。 在 他 们 的 研究 中 使 用 的 样本 都 是 距离 小 于 10 AU 和 轨道 周期 小 于 104 d 的 双星 。 
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除 此 之 外 ， 密 近 双 星 的 双星 比 还 会 随 着 有 效 温度 的 增 大 而 减 小 。Gao 等 人 "通过 分 析 SDSS 
ee ed re rr ee 
的 轨道 周期 都 在 1000 d 以 内 ，SDSS SEGUE 样本 星 的 双星 比 为 43.0% 土 2.0%，LAMOST 
样本 星 的 双星 比 为 30.0% 土 8.0%。 

Liu 研究 了 LAMOST 和 Gaia 在 太阳 邻 域 观测 到 的 恒星 的 双星 比 ， 假 设 质量 比分 布 是 
容 律 的 ， 窜 律 指数 为 Y， 发 现 双 星 比 与 7 呈现 明显 的 反 相 关 关 系 。 换 句 话 说， 具有 双星 比 的 
恒星 群 会 包含 更 多 较 大 质量 比 的 双星 。 以 y = 1.2 为 界 ， 高 y 的 恒星 质量 较 小 ， 金 属 丰 度 较 
高 ， 而 低 y 的 恒星 质量 较 高 ， 金 属 丰 度 较 低 。 高 y 的 恒星 双星 比 与 金属 丰 度 呈 反 相关 ， 与 
主 质量 则 不 相关 ， 同 时 ， 低 y 恒星 的 双星 比 与 主 质量 表现 出 明显 的 相关 性 ， 但 与 金属 丰 度 
没有 呈现 出 明显 的 相关 性 。 

主 序 星 的 双星 比 还 与 主星 质量 有 明显 的 相关 性 。Duchane 和 Kraus’ HHT EFE W 
矮星 、 主 序 前 星 和 髓 埋 的 原 恒星 的 所 有 经 验 知 识 。 对 于 太阳 型 的 恒星 和 小 质量 的 恒星 ， 轨 道 
周期 的 分 布 是 单 峰 的 。 然 而 ， 随 着 恒星 质量 的 减 小 ， 中 值 间 距 和 宽度 分 布 急 剧 减 小 ， 结 果 ， 
fE1~ 10 AU 范围 内 恒星 的 多 重 性 频率 不 会 随 着 恒星 的 质量 有 明显 的 变化 ， 双 星 比 和 其 他 
的 重要 参数 都 随 着 恒星 质量 平滑 变化 。 

鉴于 上 述 研 究 结 果 ， 我 们 利用 大 样本 的 LAMOST 数据 与 WISE 的 数据 进行 交叉 ， 解 吕 
这 一 批 双 星 的 轨道 周期 、 质 量 比 和 轨道 偏心 率 ， 并 且 分 析 了 双星 的 轨道 周期 分 布 、 质 量 比分 
布 和 轨道 偏心 率 分 布 ， 之 后 与 Raghavon 等 人 ”的 结果 进行 了 对 比 。 

本 文 结构 如 下 : 第 2 章 介 绍 求解 双星 轨道 参数 所 使 用 的 研究 方法 和 模型 ， 第 3 章 介绍 使 
的 数据 ; 第 4 章 将 此 方法 应 用 在 这 批 数据 上 所 得 出 的 结果 ; 最 后 一 章 给 出 总 结 和 展望 。 
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2 研究 方法 及 方法 验证 


2.1 ”研究 方法 

对 于 圆 形 轨道 的 特殊 情况 ， 物 体 围绕 质量 中 心 运动 的 速度 是 恒定 的 。 但 对 于 偏心 轨道 
的 一 般 情况 ， 速 度 是 位 置 和 时 间 的 连续 变化 函数 。 所 以 椭圆 轨道 中 的 速度 是 位 置 的 强 函 数 ， 
特别 是 对 于 非常 偏心 的 轨道 。 恒 星 在 椭圆 轨道 上 的 位 置 不 是 时 间 的 简单 函数 ， 对 于 椭圆 轨 
道 ， 描 述 这 种 位 置 和 时 间 的 依赖 关系 的 基本 方程 称 为 开 普 勒 方程 ”。 椭 圆 轨道 的 开 普 勒 方 
程 为 : 


BesinB= (tT) . (1) 


A (1) AY ad ee Be A a ER, Ze A TBE PE RRS EA i. EA i 
fA, e 为 轨道 偏心 率 ，P 为 轨道 周期 ，t 为 时 间 ， TT 为 积分 常量 。 并 且 ， 我 们 从 椭圆 轨道 的 
性 质 中 还 非常 容易 得 到 : 


rcosO=a(cosE—e) , (2) 


r=a(l—ecosE) , (3) 
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cos H —e 
COS 0 = I= eosi . (4) 


当 我 们 给 出 空间 中 双星 轨道 相对 于 观测 者 视线 的 方向 之 后 ， 根 据 上 述 公 式 经 过 推导 ， 非 常 容 
易 得 到 由 轨道 运动 产生 的 径 向 速度 ， 表 示 为 


2rasini 


Viad = Pe) [cos(0 +w) +ecosw] . (5) 

径 向 速度 的 最 终 表达 式 通 常 写成 ; 
ada 王 开 [cos(0 十 w) 十 ecosw| 十 yY , (6) 
2rasini 号 


其 中 K= pg aie 是 速度 曲线 的 半 振 幅 ，? 是 系统 速度 。 如 此 ， 当 我 们 知道 “组 双星 
的 参数 轨道 周期 P， 质 量 比 gs， 轨道 偏心 率 e， 轨 道 倾角 与 视线 方向 夹 角 4， 进 动 角 w, A 
位 角 0 以 及 主星 质量 m1， 就 可 以 画 出 双星 的 视 向 速度 曲线 。 

如 果 我 们 有 n 次 视 向 速度 的 测量 ， 己 知 主 质 量 m 的 情况 下 ， 我 们 用 EMCEE 软 件 包 进 
行 马尔 可 夫 链 蒙特 卡 罗 模 拟 ”。 并 且 我 们 采用 MCMC 中 随机 抽取 的 概率 最 大 值 和 15%、 
85% 百分数 作为 轨道 参数 的 最 佳 值 以 及 它 的 不 确定 性 ， 分 别 得 到 轨道 周期 Pp， 质量 比 g， 轨 
道 偏 心率 e， 轨 道 倾 角 与 视线 方向 夹 角 i， 进 动 角 w， 轨 道 零 相位 9。 
2.2 ”恒星 质量 的 确定 

我 们 通过 比较 每 颗 恒 星 的 有 效 温 度 (Tee) 表面 重力 (lg g) 和 金属 丰 度 ([Fe/H]) 与 
PAdova and TRieste 恒星 演化 模型 (PARSEC) T, jit XGboost 算法 估计 每 颗 恒 星 的 
质量 。 我 们 使 用 PARSEC 恒星 模型 提供 的 Tg, le g 和 Fe/H] 作为 训练 数据 ， 恒 星 质 量 作 
为 训练 标签 ， 然 后 训练 XGboost 模型 。 具 体 的 训练 过 程 使 用 了 均 方 根 误差 (RMSE) 作为 
损失 函数 (loss function)。 经 过 独立 的 测试 样本 进行 测试 ， 我 们 发 现 测 试 数 据 集 真实 恒星 
质量 与 预测 的 恒星 质量 之 间 的 均 方 根 误差 为 0.02 Moo 紧 接 着 ， 我 们 使 用 10 倍 交叉 验证 
(cross-validation) 确认 该 算法 带 来 的 总 体 不 确定 度 为 0.02 

此 外 ， 我 们 使 用 蒙特 卡 洛 方法 估计 了 推算 的 每 一 个 恒星 质量 的 误差 。 具 体 过 程 是 进行 
1000 次 抽样 ， 每 次 抽样 从 Tom, lg g 和 [Fe/H] 的 高 斯 分 布 中 随机 抽取 一 组 新 的 值 ， 以 这 些 参 
数 的 测量 值 作为 平均 值 ， 以 它们 的 不 确定 度 作为 标准 差 。 然 后 以 这 组 抽样 的 参数 通过 训练 
好 的 XGboost 模型 估计 恒星 质量 ， 进 行 多 次 测量 以 后 ， 我 们 采用 每 组 测量 值 的 中 位 数 作为 
推算 的 恒星 质量 ， 标 准 差 作为 恒星 质量 的 误差 。 人 恒星 质量 的 典型 不 确定 度 为 0.03 Mo (中 位 
数 )。 
2.3 ”方法 的 验证 

首先 ， 我 们 模拟 一 个 恒星 的 情况 。 假 设 一 颗 已 知 其 轨道 周期 、 质 量 比 和 轨道 偏心 率 的 假 
星 ， 然 后 估计 其 轨道 参数 ， 并 且 与 这 颗 星 预 设 的 真 值 进行 比较 ， 从 而 可 以 确定 该 方法 的 准确 
性 。 当 我 们 已 知 恒星 的 观测 时 间 和 视 向 速度 时 ， 将 随机 选择 所 有 其 他 参数 并 运行 MCMC. 
lg P, q, e 可 以 很 好 地 收敛 到 某 个 值 ( 见 图 1)， 该 假 星 的 轨道 参数 的 真实 值 为 lg P=1.90, 
q=0.7, e=0.3, mı =1M ， 观 测 次 数 为 20 次 ， 在 观测 周期 为 60 d 内 均匀 分 布 。 
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lg P=1.9010% 


q=0.69:9% 


e=0.31 +9% 


1.88 1.92 1.96 2.00 0.30 0.60 0.90 0.2 04 0.6 
lg P q e 


注 ， 这 颗 模 拟 星 的 轨道 参数 的 真实 值 为 lg P=1.90, g=0.7, e=0.3, BOA 10 次 比较 均匀 的 视 向 速度 的 测量 ， 
使 用 EMCEE 软件 包 运行 马尔 可 夫 链 蒙特 卡 罗 模拟 的 corner 图 


1 ”模拟 星 运行 MCMC 轨道 参数 的 corner 


对 于 这 样 一 颗 恒 星 ， 我 们 模拟 了 随 着 观测 次 数 的 增加 ， 轨 道 参数 的 模型 给 出 的 值 与 真实 
值 的 差别 ， 如 图 2 所 示 ， 观 测 次 数 范围 2 ~ 20 次 ， 纵 坐标 为 lg P, q, e， 黑色 虚线 是 恒星 的 
真实 值 ， 蓝 色 线 是 模型 在 不 同 观察 时 间 下 给 出 的 值 ， 红 色 是 误差 棒 。 可 以 看 出 ， 随 着 观测 次 
数 的 增加 ， 模 型 值 越 来 越 接 近 真 实 值 ， 并 且 误 差 棒 总 体 上 趋 于 减 小 。 由 于 双星 的 视 向 速度 方 
程 由 7 个 参数 控制 ， 如 果 观 测 的 数据 点 少 于 参数 的 数量 ， 将 很 难得 到 可 靠 的 拟 合 结果 ， 由 图 
中 也 可 看 出 ， 当 观测 次 数 达到 8 次 或 9 次 时 ， 模 型 值 与 真实 值 基本 保持 一 致 。 

接 下 来 ， 我 们 统计 了 50 颗 模 拟 星 的 情况 ， 结 果 如 图 3 所 示 。 这 50 颗 恒 星 的 轨道 周期 
P 的 范围 为 10 ~ 100 d, 质量 比 g 的 范围 为 0.3 ~ 0.9， 轨 道 偏心 率 e 的 范围 为 0.1 ~ 0.7, 
m =1M ， 观 测 次 数 从 4 ~ 20 次 ， 前 三 幅 图 是 模型 给 出 的 轨道 参数 P, q, e 在 不 同 观测 次 
数 下 的 误差 ， 后 三 幅 图 是 模型 给 出 的 轨道 参数 的 误差 的 均 方 根 。 可 以 看 到 ， 随 着 观测 次 数 的 
增加 ， 轨 道 参数 的 误差 在 0 上 下 浮动 ， 误 差 的 总 体 均 方 根 随 着 观测 次 数 的 增加 而 减 小 。 并 
且 ， 如 果 我 们 能 确定 轨道 周期 已 ， 将 其 作为 已 知 量 代入 模型 中 ， 其 他 轨道 参数 的 估 值 会 更 接 
近 真 实 值 。 
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注 : 这 颗 模拟 星 的 轨道 参数 的 真实 值 为 lg P=1.90, q=0.7, e 二 0.3， 观 测 次 数 从 2 次 到 20 次 ， 黑 色 虚 线 为 模 
拟 星 轨 道 参数 的 真实 值 ， 蓝 色 线 为 模型 给 出 的 值 ， 红 色 线 为 误差 棒 。 


图 2 ”模拟 星 在 不 同 观测 次 数 下 轨道 参数 模型 值 与 真实 值 的 差别 


1 8 12 16 20 4 8 12 16 20 “a 8 12 16 20 
观测 次 数 观测 次 数 观测 次 数 


5 图 3 50 颗 模拟 星 在 不 同 观测 次 数 下 模型 给 出 轨道 参数 的 误差 和 误差 的 均 方 根 


3 数 W 


LAMOST 是 一 架 新 型 的 大 视 场 兼备 大 口径 望远镜 ， 也 叫 郭 守 敬 望远镜 呈 。LAMOST 
应 用 薄 镜 面 主动 光学 加 拼接 镜面 主动 光学 技术 ， 在 曝光 1.5 h 内 可 以 观测 到 暗 达 18 mag 的 
v 天 体 ， 使 其 成 为 大 口径 兼 大 视 场 光学 望远镜 的 世界 之 最 。 同 时 ， 采 用 并 行 可 控 的 光纤 定位 技 
全 术 ， 在 5° 视 场 上 放置 4000 根 光 纤 ， 可 以 同时 获得 4000 个 天 体 的 光谱 ， 使 其 成 为 世界 上 光 
= 谱 获取 率 最 高 的 望远镜 ” 。 
C 我 们 使 用 的 是 LAMOST DR7 和 部 分 DR8 的 中 分 辨 率 数据 ， 使 用 蓝 端 的 观测 时 间 、 
观测 速度 及 速度 的 误差 ， 并 且 与 Chen 等 人 " WISE 食 双 星星 表 进 行 交叉 ， 筛 选 出 其 中 
LAMOST 光谱 的 信 噪 比 大 于 10 的 ， 观 测 次 数 大 于 8 次 的 星 ， 共 有 151 颗 。 我 们 将 这 151 颗 
星 代入 模型 中 ， 主 星 质量 m 是 由 上 述 恒星 质量 确定 的 方法 给 出 ， 将 模型 得 出 的 周期 与 
WISE 给 出 的 周期 进行 比较 。 从 图 4 可 以 看 出 ， Pwisk 值 和 ELAMOST 值 的 分 布 大 部 分 在 

Puamost PwisE 

1 左右 ， 少 部 分 为 1/2, 2 和 3， 说 明 我 们 模型 给 出 的 轨道 周期 值 与 WISE 给 出 的 基本 一 致 ， 
少 部 分 为 2 倍 周期 或 者 3 倍 周期 ， 所 以 我 们 认为 WISE 给 出 的 周期 值 是 基本 可 靠 的 。 
由 于 LAMOST 中 分 辩 率 给 出 的 视 向 速度 的 观测 点 比较 少 ， 得 到 的 周期 缺乏 可 靠 度 ， 因 
此 用 WISE 给 出 的 基于 光 变 曲线 的 周期 进行 接 下 来 的 双星 轨道 参数 的 测定 ， 我 们 使 用 WISE 
给 出 的 周期 作为 已 知 量 代 入 模型 中 ， 在 这 151 颗 星 中 挑选 出 了 视 向 速度 的 相位 分 布 比较 广 
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4 151 颗 星 模型 给 出 的 轨道 周期 与 WISE 的 比较 


1 1 mi 
6 7 8 


并 且 比 较 均匀 的 ， 速 度 误差 较 小 并 且 模 型 拟 合 出 的 速度 曲线 较 好 的 23 颗 星 。 从 LAMOST 
得 到 的 这 23 颗 星 的 有 效 温 度 (Ter). KHEN (lg g) 和 金属 丰 度 ([Fe/H]) 的 分 布 见 图 5， 


这 一 批 星 的 有 效 温 度 (Toe) 主要 集中 在 5000 ~ 7500 K 左右 ， 表 面 重力 (lg g) 主要 集中 在 


3.8 ~ 4.5 左右 ， 金 属 丰 度 ([Fe/H]) 集中 在 -0.5 ~ 0.5 左右 。 


5 6r 


4 


3 


0 


5000 6000 7000 8000 
Tor lg g 


5 ”样本 星 的 有 效 温度 、 表 面 重力 和 金属 丰 度 的 分 布 图 


4 结 R 


3.6 4.0 4.4 N95 -0.75 -0.25 0.25 


[Fe/H] 


我 们 将 LAMOST 给 出 的 23 颗 样 本 星 的 观测 时 间 与 视 向 速度 代入 模型 中 ， 同 时 将 WISE 


给 出 的 轨道 周期 已 作为 已 知 量 代 入 模型 中 ， 主 星 质量 m 依然 由 上 述 方法 确定 ， 得 到 双星 的 


轨道 参数 ， 发 现 由 于 LAMOST 的 视 向 速度 的 弥散 比较 大 ， 因 此 算出 的 轨 i 


道 参数 不 太 准 确 ， 


如 轨道 偏心 率 e 整体 较 大 ， 与 实际 不 符 。 所 以 我 们 用 章 博 等 人 改正 过 的 LAMOST 的 视 向 速 
E (已 投稿 ， 视 向 速度 是 由 交叉 星 得 到 的 )。 由 视 向 速度 曲线 进行 拟 合 得 到 了 这 23 颗 星 的 轨 
道 参数 : 轨道 周期 P、 质 量 比 g 和 轨道 偏心 率 e。 表 1 列 出 这 23 颗 样 本 星 的 轨道 参数 ， 轨 道 
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1 期 


WISEID 


表 1 23 颗 样本 星 的 轨道 参数 
IA IRA mi P q e sin i w t0 


/() /() /Mo /d /(@) /d 


J221437.8- 
J230015.44 
J222153.3- 
J040103.2- 
J204052.3+ 
J171648.3- 
J203105.2- 
J220742.2- 
J032504.44 
J064217.3- 
J061214.7- 
J171552.0- 
J081117.44 
J085709.7- 
J040934.9+ 
J153541.0- 
J112924.7- 
J090807.14 
J083747.3- 
J013437.0- 
J063546.14 
J091415.24 


+ 281723 
+ 323333 
+ 280246 
+ 202449 
+401318 
+ 544614 
+ 384659 
+ 305001 
+ 584049 
+ 201648 
+ 220119 
-535927 
-525234 
185643 
-500056 
454435 


-421220 
-110115 
-450938 
-192827 
192524 


333.657 743 28.289735 1.791 222334 0.544 663086 0.674 675 174 0.117 908281 0.467 003 711 0.385 584 144 0.307 468 756 
345.064 241 32.559297 0.964101195 0.307 457918 0.771 423 729 0.134 800612 0.661073 991 4.396 511 295 0.293 132 962 
335.47242 28.046 322 1.154322 386 0.311 734421 0.614042 531 0.178 825773 0.437 295 641 4.873 471 255 0.285 491 638 
60.263 4847 20.413 9633 1.630 368 352 0.433 540559 0.315 522 722 0.125 607 523 0.627 727929 3.285 437 746 0.237 562 117 
310.218 2339 40.221 8214 0.96293509 0.338 335 268 0.418 834 202 0.110067 878 0.980 271 389 5.365 733 657 0.190 193 223 
259.201 471 54.770 727 1.172252 417 0.357 660936 0.259 209 749 0.034 664043 0.689 417 293 3.982 060 061 0.112 464325 
307.771 8648 38.783 460 8 0.878 642917 0.291 972389 0.754171 353 0.063 771714 0.711825 28 3.743 993 253 0.127 863 504 
331.926 057 30.833772 1.168613911 0.357928 71 0.528 989 381 0.032 529128 0.687 251366 2.8726892 0.207469 489 
51.26866 58.680412 1.505527 854 1.148 404735 0.855 508 258 0.026 591 064 0.920 143 248 4.725 574 477 0.372 302 505 
100.572 397 4 20.280 096 7 1.180399656 0.38161504 0.519 390486 0.119 354769 0.771086 436 2.740 711 246 0.324 431 755 
93.061 5933 22.022 2208 1.599 668 264 0.500 346 968 0.430978 096 0.040 452 587 0.689 643011 4.109 872 996 0.251 522 679 
258.966 533 53.991001 1.018016815 0.295575 713 0.313 229 203 0.082 219 462 0.821019 562 4.858 314 482 0.104 740 702 
122.82257 52.876349 1.325948 1.426 284453 0.499 622 524 0.123 881 166 0.747563 514 6.065 476 742 0.848 111015 
134.290 445 18.9455617 1.20972681 0.28947327 0.293 928 349 0.077 281358 0.682 556 723 3.248 759 612 0.133671 144 
62.395 4698 50.015 8429 1.851095 796 0.533 755 852 0.478 911 282 0.080 740 754 0.894907 476 4.681 461 401 0.146 157 355 
233.920 788 9 45.743 2648 1.155 333996 0.337 701936 0.156 531 208 0.073 787 444 0.104 889 245 2.964 687 634 0.240 060 529 


+010046 172.353 054v8 1.0129731 1.368018 389 0.644 216v466 0.427 316 985 0.111657 633 0.492 572089 3.570 923 478 0.290 984 949 


137.03 42.20559 1.174523711 0.460 734182 0.719 662 643 0.096 506v832 0.908 854 483 5.824 416 701 0.272 872 997 
129.44717 11.020942 1.161530495 0.358 980788 0.498 87964 0.043 393 804 0.749 191575 1.741 300 403 0.223 688 875 
23.654478 45.160768 1.041 315v079 0.434 317927 0.344549 981 0.034879 339 0.879 507 178 0.252 046 384 0.263 450 885 

98.942 5608 19.474457 1.274019241 0.479 583034 0.624 909 935 0.065 323 159 -0.296 568 863 2.412 327 744 0.304 982 161 
138.563 3736 19.423 5692 1.17110002 0.418 881945 0.196 743 916 0.089 797 547 0.745 717 903 0.766 901 326 0.297 817 684 
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周期 由 WISE 基于 光 变 曲线 的 周期 得 到 ， 质 量 比 q 如 表 1 所 示 ， 并 且 它 的 平均 不 确定 度 约 
为 0.1， 轨 道 偏心 率 e 的 平均 不 确定 度 约 为 0.02。 图 6 展示 了 这 23 颗 星 代入 模型 中 得 到 的 轨 
道 参数 对 视 向 速度 曲线 的 拟 合 情 况 。 轨 道 周 期 P、 质 量 比 g 和 轨道 偏心 率 e 的 分 布 情况 如 
图 7 所 示 ， 可 以 看 出 ， 轨 道 周 期 已 是 由 WISE 给 出 的 ， 大 部 分 小 于 1 qd， 都 是 短 周期 的 密 近 
双星 ， 分 布 呈 现 出 肾 律 分 布 的 趋势 ， 质 量 比 的 分 布 整体 比较 平 ， 轨 道 偏心 率 e 由 于 动力 学 测 
汐 作 用 造成 的 轨道 贺 化 ， 所 以 短 周期 的 密 近 双星 的 轨道 偏心 率 都 比较 小 ， 模 型 给 出 的 轨道 偏 


心率 也 如 预期 ， 整 体 偏 小 ， 都 在 0.2 以 内 。 


Yaa/ km's ) 


6 23 颗 样 本 星 的 视 向 速度 曲线 拟 合 结果 


接 下 来 ， 将 我 们 得 出 的 结果 与 Reghavon 等 人 ”的 结果 进行 对 比 。 我 们 画 出 了 轨道 周期 


与 质量 比 及 轨道 偏心 率 之 间 的 关系 图 (如 图 


8 所 示 )， 横 坐标 是 轨道 周期 取 10 的 对 数 lg P, 
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7 ”样本 星 的 轨道 周期 、 质 量 比 和 轨道 偏心 率 的 分 布 图 
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纵 坐 标 分 别 为 质量 比 g 和 轨道 偏心 率 e， 在 轨道 周期 lg P 为 一 0.6 ~ 0.2 范围 内 ， 质 量 比 g 
FEA tie SES 0.1 ~ 0.9 的 空间 ， 并 且 整 体 趋向 于 比较 大 的 值 ， 但 是 轨道 偏心 率 e 都 小 于 


0.2， 并 且 大 部 分 集中 在 0.1 左右 。 
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图 8 ”轨道 周期 与 质量 比 及 轨道 偏心 率 的 关系 


之 后 我 们 将 图 8 的 轨道 周期 与 质量 比 和 轨道 偏心 率 分 别 与 Reghavon 等 人 ”的 结果 进行 
比 的 关系 ， 空 心 圆 是 我 们 的 23 颗 短 周期 样本 星 ， 加 号 、 


对 比 。 图 9 给 出 ， 轨 道 周 期 与 质量 


实心 三 角形 和 空心 正方 形 分 别 表示 Reghavon 等 人 村 


本 中 的 双星 系统 、 三 星系 统 和 四 星系 


统 ， 可 以 看 出 我 们 的 结果 与 Reghavon 等 人 中 的 结果 基本 一 致 ， 质 量 比 基 本 覆盖 全 空间 。 


图 10 给 出 轨道 周期 与 轨道 偏心 率 的 关系 ， 同 样 ， 


空心 圆 是 我 们 的 23 颗 短 周期 样本 星 ， 


加 号 、 实 心 三 角形 和 空心 正方 形 分 别 表示 Reghavon 等 人 ”样本 中 的 双星 系统 、 三 星系 统 和 
四 星系 统 。 曲 线 虚 线 代 表 一 个 偏心 率 的 边界 ， 在 其 左边 ， 轨 道 过 于 扇 ， 一 颗 伴星 到 另 一 颗 伴 
星 近 星 点 的 时 候 在 1.5Ro 之 内 ， 因 此 很 可 能 发 生 碰撞 。 我 们 的 样本 星 的 轨道 偏心 率 都 很 小 ， 


都 小 于 0.2， 在 曲线 内 。 
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我 们 的 23 颗 短 周期 样本 星 ， 加 号 、 实 心 三 角形 和 空心 正方 形 分 别 表示 Reghavon 等 人 四 样本 中 的 
双星 系统 、 三 星系 统 和 四 星系 统 。 


注 ， 空心 
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图 9 ”轨道 周期 与 质量 比 的 关系 
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统 、 三 星系 统 和 四 星系 统 。 
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5 总 结 与 展望 


我 们 使 用 了 LAMOST 中 分 辨 率 的 DR7 和 部 分 DRS 的 数据 ， 得 到 了 23 颗 样 本 星 。 通 
过 其 视 向 速度 测定 了 其 轨道 参数 ， 轨 道 参数 基本 符合 预期 ， 结 果 如 表 1 所 示 。 除 此 之 外 ， 还 
存在 一 些 问题 。 对 于 我 们 轨道 参数 的 结果 ， 尽 管 轨 道 偏心 率 e 都 在 合理 的 范围 内 ， 但 还 是 偏 
大 一 些 ， 特 别 是 图 6 中 第 一 列 的 倒数 第 二 张 图 ， 轨 道 偏心 率 还 是 比较 大 。 这 可 能 是 由 视 向 速 
度 测量 的 误差 造成 的 或 者 观测 时 间 的 采样 相位 不 均匀 造成 的 ， 这 类 问题 希望 后 续 通过 对 数 
据 处 理 的 精度 提高 和 观测 次 数 的 增加 来 改善 。 未 来 ， 这 类 方法 经 过 扩展 可 以 变 成 一 个 层次 化 
的 贝 叶 斯 模型 ， 在 不 做 单 颗 星 的 轨道 参数 的 条 件 下 还 可 以 估计 轨道 的 性 质 。 
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Determination of Orbital Parameters of Binary Stars Based 


on Spectrum and Photometric Data 
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Academy of Sciences, Beijing 100049, China) 


Abstract: Binary star systems or multi-star systems are a common result of stellar evolu- 
tion. Therefore, the dependence of the stellar multiplicity and other properties of the binary 
star system on the primary mass and the stellar formation environment play important role 
in stellar evolution. There are many ways to study binary stars, such as the light curves 
method, the velocity curves method, and the Hertzsprung-Russell diagram method. In re- 
cent years, with the emergence of many large-scale astronomical survey observations, we 
have obtained a large amount of photometric and spectral survey data. For these large-scale 
survey data, many statisical methods have been developed to study the properties of binary 
stars. We use the WISE eclipse binaries with the medium-resolution spectra of LAMOST 
DR7 and part of DR8, and obtained 151 short period solar-type stars, and then selected 
those with more than 9 observations and relatively uniform phase distribution. By ana- 
lyzing their radial velocities, the orbital parameters such as the orbital period, mass ratio, 
orbital eccentricity of this sample are obtained. We find that the period we obtained consis- 
tent with that given by WISE. Thus, we use the period of WISE as a known quantity to solve 
for other orbital parameters, which makes other orbital parameters more accurate. Among 
them, 23 stars were selected and the radial velocity curves are well fitted by binary models. 
The statistical properties of these parameters are consistent with the results of Ragavan et 
al. For stars with orbital period of less than 10 days, their mass-ratio basically covers the 
entire parameter space from 0.2 to 1, while their orbital eccentricities are very small, all 


smaller than 0.2. 
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